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Universality of the saturation scale and the initial eccentricity in heavy ion collisions 
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Recent estimates that Color Glass Condensate initial conditions may generate a larger initial 
eccentricity for noncentral relativistic heavy ion collisions (relative to the initial eccentricity assumed 
in earlier hydrodynamic calculations) have raised the possibility of a higher bound on the viscosity of 
the Quark Gluon Plasma. We show that this large initial eccentricity results in part from a definition 
of the saturation scale as proportional to the number of nucleons participating in the collision. 
A saturation scale proportional to the nuclear thickness function (and therefore independent of 
the probe) leads to a smaller eccentricity, albeit still larger than the value used in hydrodynamic 
models. Our results suggest that the early elliptic flow in heavy ion collisions (unlike multiplicity 
distributions) is sensitive to the universality of the saturation scale in high energy QCD. 

PACS numbers: 24.85.+p, 25.75.-q, 12.38.Mh 



I. INTRODUCTION 

Data from RHIC on bulk properties 0, S H 3 of the 
matter produced in heavy ion coUisions, is quahtatively 
described by hydrodynamic calculations . Increasingly 
precise experimental data now enables us to probe details 
of initial conditions for the hydrodynamical calculations. 
These initial conditions may be obtained from the Color 
Glass Condensate (CGC) framework 6] describing the 
phenomenon of parton saturation (2lil2i in high energy 
nuclear wavefunction. 

One particular application of the CGC framework to 
the early stages of heavy ion collisions has been the for- 
mulation of the problem as a solution of classical field 
equations 0, 0, lisll a nd the numerical solu- 

tion of these equations01ll El El, 111 III 13, El 
I2I I2I 0) H^- By solving classical equations, one ef- 
ficiently sums all leading order Feynman graphs con- 
tributing to the inclusive multiplicity. We shall hence- 
forth refer to this approach as the Classical Yang-Mills 
(CYM) approach. An alternative approach is to assume 
the framework of k± factorization, where the inclusive 
multiplicity of produced gluons is expressed as a convo- 
lution of unintegrated gluon distributions from the pro- 
jectile and the target. While exact to leading order for 
proton- nucleus collisions [2^ l28l |. this approach is 
only approximate in heavy ion collisions at central rapidi- 
ties. A popular version of the k± factorization approach 
applied to heavy ion collisions, with a particular ansatz 
for the unintegrated gluon distributions, is the KLN ap- 
proach m Ha m lH 113. The solution of the CYM 
equations incorporates, to lowest order in the coupling, 
all terms that violate k± factorization; these solutions 
however are numerically intensive. On the other hand, 
the KLN ansatz, while approximate, provides an analytic 
expression that captures some key features of saturation. 

A striking signal of collective hydrodynamical behav- 
ior of the matter produced at RHIC is elliptic flow 
I35, 36, 37] . The large elliptic flow observed at RHIC re- 
quires strong interactions during the first fermis of the 
collision. The elliptic flow in this strongly interacting 



system is particularly sensitive to the initial eccentricity 
(the anisotropy of the energy density in the transverse 
plane) , 



(1) 



where e is the local energy density ^. In a hydrodynam- 
ical description of heavy ion collisions, the eccentricity 
of the initial condition is the primary factor determining 
the elliptic flow observed in the final state [39l| . 

In this short note, we will study how the eccentricity 
is computed in the CGC framework in the CYM and 
KLN models. We shall first discuss different ways used 
in the literature to understand the transverse coordinate 
dependence of the saturation scale, in particular based 
on its universality. We shall then explain how these lead 
to different results for e, in particular the recent results 
in Refs. 0,1^. (See also Ref. 42] for a discussion of 
related issues for uranium- uranium collisions.) 



II. A UNIVERSAL SATURATION SCALE 

In the CGC framework, the production of particles 
in the initial stage of a heavy ion collision is controlled 
by one parameter, the saturation scale Qs- The initial 
spatial anisotropy of the energy density must thus come 
from the transverse coordinate dependence of the satu- 
ration scale. In the original MV model [isL IS lisf the 
large x degrees of freedom are thought of as indepen- 
dent static (independent of light cone time) color charges 
with a transverse color charge density g^^J? ■ The den- 
sity of color charges is a property of the nuclear wave- 



^ Alternately, the eccentricity i s de fined as the anisotropy of the 
entropy density in some works l38l . As we will see, this definition 
is less sensitive to different CGC initial conditions. It is however 
not the definition used in comparisons of detailed hydrodynamic 
models to RHIC data. 
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function and independent of the probe — it is univer- 
sal If the color charge density arises as a superpo- 
sition of independent large x partons, it is natural that 
the density should be proportional to the nuclear thick- 
ness function Ta, thereby providing an average measure 
of the number of these fast partons. It was later shown 
m EUS m El El that the MV model does in- 
deed exhibit saturation, with Qs ~ g^n, and thus also 
the saturation scale Qs should be universal with 



A 



■7a (xt) 



(2) 



Here Qso is the average saturation scale, Ta(xt) — 
dz p{r) is the thickness function, xt is the coor- 
dinate relative to the center of the nucleus and p(r) is 
the Woods-Saxon nuclear density profile, normalized as 
/ d'^rp(r) = A. It has also been argued [H3 that at 
sufficiently high energy the dependence of the saturation 
scale on Ta and thus A would change due to high energy 
evolution, but we will not consider this possibility here. 

The classical field picture of the MV model, with this 
universally defined saturation scale, was applied in both 
analytical iS. .10.^ 11 .. 1 2. 13) and numerical classical Yang- 

Mills (CYM) TaMMiBlilliijpiIiiiia 

computations of the "Glasma" '54!| fields in the initial 
stages of the collision, Nonccntral coUisions^fith finite 
impact parameter were studied in Refs. |l 9l[2l | and par- 
ticularly in the detailed study of Ref. l2d . 

In the "KLN" approach, H ES IH El El, par- 
ticle production is computed using unintegrated gluon 
distribution functions depending on the saturation scale 
Qaix, xt). These calculations used a saturation scale de- 
pendent on the number of participant nuclei: 



■(xt) ~iVpl,t(xT), 



(3) 



with 



7Vpt,(xT)=T^(xT) 
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Note that as this form involves the thickness functions of 
both nuclei A and B, it is not manifestly universal. While 
a nonuniversal saturation scale has been used in many 
phenomenological studies, such as the EKRT "final state 
saturation" model of Ref. [5|, the logic for introducing 
the second term in the original KLN works was not the 
final state picture. The A^pait definition was used instead 
to model the interaction probability for the scattering 
in a manner as close to the original Glauber framework 
as possible. There arises however the danger of double 
counting the interaction probability because this proba- 
bility is already taken into account in the derivation of 
k± factorization ^6, 27, 28]. The dan gers of overcount- 
ing are avoided in the universal prescription for Qg in 
Eq. 




FIG. 1: The area in the transverse plain that is most impor- 
tant for our discussion here is marked with the gray oval, near 
the edge of nucleus A and the center of nucleus B. The same 
discussion applies for the edge of nucleus B near the center 
of nucleus A. 



We must emphasize that the choice between using the 
k± factorized formalism and solving the classical field 
equations does not dictate how the saturation scale de- 
pends on the transverse coordinate. There is no technical 
impediment to using (Q^)^ ~ Ta in the k± factorized 
formulation or (Q^)^ ~ ^pmt ™ f^e classical field com- 
putation. In the numerical study using the classical field 
equations in Ref. both prescriptions, Eqs. Q and Q 
were considered and were found to give similar results 
for the impact parameter dependence of the multiplicity. 
However, as we shall show in the following, they do yield 
different results for the eccentricity. 



III. HOW THE CHOICE OF INFLUENCES 
THE ECCENTRICITY 



For computing the eccentricity, the crucial contribu- 
tions (in the geometry of the transverse plane of the scat- 
tering) come from regions where the saturation scale in 
one nucleus is significantly larger than in the other. Let 
us first review what the spectrum of gluons looks like in 
^4) this case, with two saturation scales Qsi and satis- 
fying Qsi < Qs2 iil2i>i28j. Parametrically, the spectrum 
of produced gluons behaves as 
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Integrated over transverse momenta, this gives 
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FIG. 2: The eccentricity as a function of impact parameter. 
The classical field CGC result with two different infrared cut- 
offs m is denoted by CYM. The traditional initial eccentric- 
ity used in hydrodynamics is a linear combination of mostly 
"Glauber A^part" and a small amount of "Glauber A'^coii". The 
"KLN" curve is the eccentricity obtained from the CGC cal- 
culation in Refs. |40 l l4l | . 



neglecting logarithmic corrections ^ In ((5s2/*3si)- The 
additional dependence on Qs2 in the transverse energy, 
relative to the multiplicity, holds the key to the following 
discussion. 

The difference between the two definitions of the trans- 
verse coordinate dependence of the saturation scale, 
Eqs. |(2J) and ijSJl is the largest in the region near the 
edge of one nucleus (labeled as nucleus A) and in the 
center of the other (nucleus B), so that < Qf; the 
geometry is illustrated in Fig. ^ The smaller saturation 
scale approaches zero as (Q^) ~ Ta regardless of the 
definition of Qs (Eq. (|2J) or Eq. Q). But the behavior 
of the larger saturation scale Qf is different in the two 
cases. Using the universal definition of Qs in Eq. ^ 
Qf is large, (Qf)^ ~ Tb- In contrast, the non-universal 
.^part-definition of Qs in Eq. (j2Jl suggests that the larger 
saturation scale also approaches zero as aMNTATs- 

Because the multiplicity © only depends on the 
smaller saturation scale Qf, the difference in the gluon 
multiplicities between the two definitions Eqs. Q and Q 
is small. This explains the numerical observation in 
Ref. [231 that both the KLN prescription for Qs and the 
universal CYM one give very similar results for the cen- 
trality dependence of the multiplicity. The larger sat- 
uration scale Qf and therefore the energy density are, 
however, very different in the two cases. This difference 
is accentuated in the eccentricity . With the A'^part def- 
inition (O, the energy density in this edge region is sup- 
pressed relative to the universal definition in 10), thereby 
leading to a larger eccentricity. 

The eccentricities obtained using the different trans- 
verse coordinate dependences of the saturation scales are 



shown in Fig. |5] The CYM eccentricity in the plot is 
calculated a,t t — 0.25 fm, while the KLN result does not 
depend on time. The KLN A^part definition of Qs leads 
to the largest eccentricity. The universal CYM definition 
gives smaller values of e albeit larger than the traditional 
parametrization (used in hydrodynamical model compu- 
tations) where the energy density is taken to be propor- 
tional to the number of participating nucleons. This re- 
sult is also shown to be insensitive to two different choices 
of the infrared scale m which regulates the spatial ex- 
tent of the Coulomb tails of the gluon distribution. We 
observe that the values of e from the CYM computa- 
tion are close to those obtained from an energy density 
parametrization following binary collisional (A^coii) scal- 
ing. This result can be explained qualitatively as fol- 
lows. In the classical Yang-Mills calculation the total 
multiplicity of gluons scales as Qg, where Qs is the domi- 
nant transverse momentum scale of the produced gluons, 
depending on both saturation scales Q^ and Qf . The 
multiplicity of produced gluons ~ Qg turns out to be 
roughly proportional to iVpart- The energy density, on 
the other hand, scales as Qs, and one expects it to scale 
as (A'^part)''' with some 7 > 1. It is therefore natural to 
expect the eccentricity in a saturation model to be larger 
than the traditional one following from A/'part-scaling of 
the energy density. However, we see no reason in general 
for it to exactly mimic the result from A'coii-scaling. 

In Fig. O we show a plot of the centrality dependence 
of the multiplicity for g^/i =1.6 GeV corresponding to an 
estimated gluon multiplicity of ^ 1000 in central Au-Au 
collisions at RHIC ^. The universal Qg ~ Ta prescrip- 
tion captures the observed centrality dependence of the 
multiplicity distribution. It has been argued [57l | that 
in a realistic Monte Carlo implementation the KLN for- 
malism can be recast in a form where the multiplicity is 
equivalent to one calculated from universal unintegrated 
gluon distributions. It appears unlikely however that this 
equivalence holds for other observables. 



IV. CONCLUSIONS 

We have shown in this brief note that the initial eccen- 
tricity of a relativistic heavy ion collision, computed in 
the Color Glass Condensate framework, is very sensitive 
to the transverse coordinate dependence of the satura- 
tion scale Qs. When Q^ is proportional to the number of 
participants A^part i the energy density produced (near the 



^ The valu e g^^ = 1.6 differs from the previous estimate of 2 
GeV [ifij l primarily because these estimates had very low infrared 
cut-offs of order m ~ 1/ Rj^ . For finite nuclei an infrared scale m 
of the order of the surface diffuseness of the Woods-Saxon den- 
sity profile is required to regulate the Coulomb tails of the gluon 
field at large distances. While the dependence on m is weak, 
changing it by a factor of 10 does change the best estimate for 
by ~ 20%. 
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FIG. 3: The centrality dependence of the multiphcity. The 
closed circles are the PHOBOS data for the charged particle 
mulciplicity per participant pair from Ref. [H^. The open 
squares are the CYM result for gluon multiplicity for g^fj, — 
1.6 GeV and m = 0.5 GeV and scaled by TVch/A^tot ~ 2/3. 



edge of one nucleus while near the center of the other) 
is small, leading to a large eccentricity. An argument 
based on the universality of the saturation scale leads to 
in nucleus A being proportional to the nuclear thick- 
ness function Ta of that nucleus alone. In this latter 
case, the saturation scale in the center of one nucleus 
does not depend on the other. While the two definitions 
have only logarithmic differences for multiplicity distri- 
butions, the energy density is larger and the eccentricity 
smaller in the universally scaling form of Qs than for 
the non-universal case. In both cases, the result for the 
eccentricity is larger than the values typically employed 
in hydrodynamical models which assume ideal hydrody- 
namics. Conceptually, the universality of the saturation 



scale is important because this universality is essential for 
a reliable framework to compute the properties of QCD 
at high energies. 

The larger eccentricities from the CGC initial condi- 
tions are of great phenomenological interest because they 
allow for the possibility that there may be significant vis- 
cous effects in the quark gluon plasma. Our discussion 
should therefore be taken into account in future attempts 
to place an upper bound on the viscosity of the plasma. 
The reader should note however that these eccentricities 
were computed at very early times (r « 0.25 fm) after 
the collision. After this time the eccentricity will imme- 
diately start to decrease because the edges of the sys- 
tem are expanding into the surrounding vacuum. Even if 
thermalization were to occur at the extremely early times 
of 0.6 fm (as some hydrodynamic models assume), the 
eccentricity will have decreased from the CYM values al- 
ready at 0.25 fm 0.6 fm, which is the time at which one 
needs the initial conditions. Therefore the initial eccen- 
tricities being assumed currently for hydrodynamic sim- 
ulations may not too implausible leaving open the possi- 
bility that one indeed might have a "perfect fluid" |^ . A 
fuller understanding (and therefore a firmer bound on the 
viscosity) will require a dynamical understanding of how 
eccentricity decreases from the initial "Glasma" stage to 
the thermalized stage of the heavy ion collision. 
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